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Herbert Shaw, Royal College of Science, South Kensington, S.W. 
(proposed by A. Fowler); and 

James Henry Worthington, Student in the University of Oxforc}, 
Bindon, Wellington, Somerset (proposed by H. H. Turner). 


Sixty presents were announced as having been received since the 
last meeting, including, amongst others:— 

E. T. Whittaker, The Theory of Optical Instruments, presented 
by the author; framed photograph of 14th century clock restored 
by J. J. Hall, presented by Mr. Hall; Calendrier Perp^tuel, pre¬ 
sented by M. Bosson. 

^ Seven transparencies of Sun-spots, Comet, Planets, and Nebulae, 
from negatives taken at the Royal Observatory, Greenwich, pre¬ 
sented by the Astronomer Royal. 

J. 0 . Kapteyn, On the number of stars of determined magni- 
* tude and determined galactic latitude (Pub. Astron. Laboratory, 
Groningen, No. 18), presented by Professor Kapteyn. 


On the Lunar Inequalities due to Planetary Action. 

By Ernest W. Brown, D.Sc., F.R.S. 

1. I have lately finished the computation of the terms in the 
Moon’s motion due to the actions—direct and indirect—of the 
planets, and it is of interest to compare the results with those of 
previous investigators. There are two extensive lists in print: 
one given hy Radau* in 1892, and the other by Newcomb t 
a few months ago. Both of these give only the terms in longitude, 
and I shall ^therefore limit the results set forth here for discussion 
to these terms, reserving the complete list for a memoir, containing 
also my methods and the details, which I hope to present to the 
Society at a later date. I have also omitted here the secular terms 
and those which depend solely on the lunar arguments, since the 
results for them are not in doubt, and since they will ultimately be 
combined with terms arising from the figure of the Earth. Further, 
I only consider quantities of the first order relative to the planetary 
masses, while Newcomh has included the mutual perturbations 
in his values for the direct action. But the latter do not affect 
any of the terms except one or two of long period, mentioned in 
No. 10 below. 

The great majority (about 400) of the terms are of short period, 

* “ Recherches concernant les In4galittis planetaires du mouvement de la 
Lune,” Paris Obs. Ann. (J 16m.), vol. xxi; 

t ‘ ‘ Investigation of Inequalities in the motion of the Moon produced by 
the action of the planets,” Carnegie Inst. Publ., 72. 
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and are most conveniently added to the true longitude. The terms 
of very long period are best added to the mean longitude, and 
they are so set down here • but the changes in the other elements 
for these long period terms are not given, since they can be ex¬ 
hibited as simple ratios to the mean longitude portions. Radau 
adds these long period terms to the true longitude, and I have 
therefore dropped the corresponding short period terms which he 
gives. * 

In his final results, Newcomb drops coefficients less than 
o"*oo3. I have done the same in transcribing from my manu¬ 
scripts, except occasionally where a smaller coefficient comes in the 
middle of a long series, or where Newcomb or Radau give a 
coefficient greater than o ,,- oo2. All my coefficients have been 
computed to o"’ooi except those of a few long period terms in ; 
these exceptions are noted by a comma placed before the zero, 
which is set down as the last figure in order not to break the 
continuity of the representation. 

2. The notation for the arguments of the tables is as follows :— 

l , 2D, F, the mean anomaly, the argument of the 
“variation” and the principal argument 
of the latitude of the Moon. 

V *, the Moon’s true longitude. 
w x , the Moon’s mean longitude. 

Q, V, T, M, J, S, the mean longitudes of Mercury, Venus, the 

Earth, Mars, Jupiter, and Saturn. 
h = D — F + T, the mean longitude of the Moon’s node. 

The results are arranged in six columns. The first pair, headed 
B, contains the coefficient and angle as computed by myself; the 
second pair, headed N, those of Newcomb; the third pair, headed 
R, those of Radau. The coefficients are given in units of o"'ooi. 
The angle has, except in Swj, been divided into four parts for 
convenience : <j> contains the lunar arguments l , D, F, h; a, the 
constant part of the angle so taken that the coefficients are all 
positive; j,f the multiples of T-P and T (or j , j" the multiples 
of P — T and P) so that / (or j") represents the order of the term 
with reference to the eccentricities and inclination of the Earth and 
planets. The values off (or j") and <f>, once given, are understood 
to be the same for every following term until new values are set 
down. The differences between the results of Newcomb and 
myself are marked with the letters E x ,' E 2 , E 3 , U whenever they 
exceed o" , o2o, and they are considered in No. 8 below. 

* An italic capital to distinguish from the mean longitude of Venus. 
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3. In order that the discussion of the various differences may be 
made clear, it is necessary to state certain facts in the theory of 
these inequalities which have considerable effect on the final results, 
and on the methods of obtaining them. 

The disturbing function R is of the form— 

R = S{A 0 cos Q + Aj cos (l + Q) + A 2 cos (2D -1 + Q) 

+ A 8 cos (2D +1 + Q) + A 4 cos (2D + Q) 
+ Ag cos (2D — 2I + Q) + ....}, 

where Q here and throughout this discussion is an angle depending 
only on the solar and planetary arguments, and where the summa¬ 
tion sign refers to the different arguments Q. 

Let w x be the Moon’s mean longitude, and c any one of the other 
five lunar elements: When the value of R has been substituted in 
the equations of variations and the latter solved we obtain— 

8^ 1 = S{B 0 sin Q +Bj sin(Z + Q)+.}, 

Sc = 3 { C 0 sin Q + Cj sin (Z + Q) + .... 

l cos cos J 

the term in Sc being a sine or cosine according as c is or is not an 
angular element. 

The values of the variations of the elements have to be 
substituted in— 

8F= 8 iv x + - 1Vwq + ^Sc . 

\dw 1 ) de 

• Now F consists of a non-periodic term w x and periodic terms, 
say, 

F=w 1 + w l 8inZ + v 2 sin (2D - Z) + v 3 sin (2D + Z) 4- . . . ., 

and from this we obtain easily the values of the deviatives of F 
The terms which arise from the substitution of the variations of the 
elements in 8F may therefore be divided into two classes. The 
first of these, which I have called elsewhere* the primary terms, 
consists of those terms in 8 F which arise from the variation of the 
non-periodic term; the secondary terms are those arising from the 
substitution in the remaining, that is, the periodic terms of V. Thus 
the values of Swq constitute the primary terms, and all the other 
portions of S F the secondary terms. 

Now the primary terms have the same arguments as the terms 
in R, that is, the arguments 

Q, Z + Q, 2D - Z + Q, . . . ., 
while the secondary terms have the arguments 

Q ± (Z), Q ± (2D — Z), Q ± (2D + Z), .... 

Z + Q i (Z), Z + Q + (2D — 1 ), Z + Q + (2D + Z), . . . ., 


* The Inequalities in the motion of the moon due to the direct action of the 
planets, Pitt Press, Cambridge, 1907.. 
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the arguments arising from the derivatives of V being enclosed in 
brackets. It is therefore obvious that the arguments Q * l + Q, 
. , . . which are present in the primary terms will also appear 
amongst the secondary terms; they, are due to the separation of 
a product of a sine and a cosine into the sum or difference of 
two sines. 

4. The first result to be noted is the fact that the 'sum of the 
secondaries with arguments l + Q — (l) , 2D - l + Q - (2D — l) , . . . . 
that is, those having the argument Q, is very small compared xoith 
the coefficient of the primary with argument Q.* In many cases 
these secondaries are of the same order of magnitude as the 
primary of the same argument, the chief of them in general being 
those arising from the two arguments just written down. This 
theorem is a consequence of the method of the variation of arbitrary 
constants. If we had adopted the straightforward method of 
integrating the original equations for the Moon’s motion with the 
additional value of ft due to planetary action, the terms inde¬ 
pendent of l in the coordinates arising from the terms in ft 
which contain l, and which therefore have the factor e (the lunar 
eccentricity), would have possessed the factor e 2 , while the principal 
terms due to R = A 0 cos Q would not have this factor; the former 
terms must therefore be quite small. 

This fact constitutes a useful test of the general accuracy of 
a large portion of the work, and it was satisfied in all cases as 
closely as could be expected.. It is illustrated in the following 
table for a few terms of the indirect action of Venus in longitude. 
The first column gives the arguments of the primaries (that is, of 
the terms in R), and in line with them are the resulting secondaries 
for the values of <j> which stand at the head of each column. The 
last line but one gives the sums of these secondaries, and the last 
line the coefficients of the primaries. The sums of these two lines 
constitute the complete coefficients for the indirect action. The 
constant angle to be added to each argument is omitted, as it is 
,(in these and in most cases) very nearly the same for each of the 
numbers in a given column. 


8F=o" , ooiCsin(<£ + a),<£ = fT — i'V, Venus. Values of C. 


Arg. of 
primary. 

T- 

-v 

T-2V 

2T- 

- 2 V 

3 T- 

■ 2 V 

3 T- BY 

4 T- 

- 3 V 

5 T- 

■ 3 V 

0 +1 

- 

47 

- 7 

- 

92 

+ 

18 

-13 

+ 

21 

+ 

5 

<p — l 

- 

8 

... 

+ 

21 

+ 

7 

+ 4 

- 

3 


7 

< f > + 2D — l 


42 

+ 6 

+ 

70 
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20 

+ 9 
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l8 
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<(> - 2D + l 

+ 

18 


+ 
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... 
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3 
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7 

< j > + 2 l - 2D 
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... 


... 



+ 1 





< f >- 2 l + 2D 
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2 

... 
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2 








Sums 
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8 

- 1 

+ 

4 

- 

3 

+ 1 

- 

3 


0 

<t> 


35 2 

-14 

- 

507 

+275 

-51 

+ : 

155 

+ 129 


* This result was not noticed until the computations had been practically 
completed, when the approximate vanishing of the sums of the secondaries 
in all oases pointed to a general theorem. 
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5. At first sight it would appear that the argumerits of the two 
terms formed, for example, from sums of products like siu (Z + Q) 
cos( 2 D- 7 ), that is, Z + Q ± (2D — /), would have coefficients in 
SF of the same order of magnitude. It is not so whenever the 
arguments in R and F contain Z, unless the argument in R is of 
very long period (greater than a score or so of years). This fact 
might have been predicted from considerations of the same, nature 
as those in the previous paragraph. For example, an argument 
2Z.+ Q 4- (Z) = Q + 3Z would not be expected from a direct method 
of treatment to have a coefficient nearly so large as that of 2Z + Q 
- (Z) = Q + Z. The latter is frequently sensible, while there are 
no terms with the former argument. 

This arises, again, from the peculiarities, of the method of the 
variation of arbitrary constants. In the volume referred to above 
I have shown directly from the method that if a term of R contains 
il in its argument (Z= 1, 2,3,...) and V a term i'l(i' — 1, 2, 3,...), 
then the secondary arising from the sum of these arguments is 
always very small compared with that arising from the difference, 
unless the period of the primary is very long. This result was 
proved by inserting the terms due to SI, Se, as given by the 
equations of variations, in SF. It enabled me to abbreviate the 
computations by some twenty or more per cent., since it rendered 
the separate computation of Se unnecessary except in a few easily 
recognised cases. 

6. One other fact should be noted., Although the coefficient 
of the evection (arg. 2D —Z) is of order m( 1/13) compared with 
that of the principal elliptic term (arg. Z), the numerical factors 
which it contains are so large that in certain functions of the 
lunar co-ordinates it becomes equal to or of more importance than 
the principal elliptic term. It is therefore never safe to neglect 
the former terms unless we know that the latter are quite insensible; 
This shows most strongly in the arguments Z + <£ in R for the 
indirect action. The terms which multiply SF' (the portion of 
the Earth’s longitude due to the action of the planets) are chiefly 
those due to the evection, while those which multiply SF (the 
planetary portion of the Earth’s radius vector) are due to the evec¬ 
tion and the principal elliptic term. Yet the former portions in 
several cases are larger than the latter. 

7. The differences B — R. Each of Radau’s coefficients depend¬ 
ing on the action of Yenus is to be diminished by 1/51 of its 
amount on account of the difference between the value for the mass 
of Yenus adopted by him and that adopted by Newcomb and 
myself. Each coefficient due to Mercury is to be diminished by 
7/60 of its amount for my adopted value of the mass of that planet. 

The equations of variations used by Radau appear to require 
a correction which would increase each coefficient in the ratio 
96 : 95, as I have shown in the volume previously referred to; 
This is confirmed by Newcomb’s results, for he obtains* values for 
the principal terms the same as those which I found in the volume 
quoted. : :■ 
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Iu the primary terms with arguments independent of the lunar 
angles, Radau omits the portions due to the term of R which 
depends on r 2 cos 2(V — T). These portions constitute about one- 
eighth ofj the principal portions for the indirect action, and his 
coefficients for this portion should therefore, be diminished in the 
ratio 7: 8. . The corresponding correction for the direct action is 
smaller, and its amount depends on the particular term under coht 
sideration. 

In the short period terms of the disturbing function for the in¬ 
direct action of the form l + <j> he omits the portions multiplied by 
SF'. These, as .stated in No. 6 above, constitute a considerable 
fraction of the whole. The terms proportionately most affected are 
those with this argument, but the differences for the terms of 
argument l + <f> — ( 1 ) are much larger, for the latter (secondary) terms 
are about three times the former (primary) terms. 

. The discovery of the causes of divergence between the results of 
Radau and myself was rendered difficult owing to the fact that, 
although his secondary inequalities for the Venus terms of argument 
<fi are given separately, and, are in several cases quite large when 
the theorem of No. 4 above shows that they ought to be quite small, 
yet his final coefficients vdry nearly agreed with mine. As a matter 
of fact, I have found, from an examination of the .separate portions 
of my results, that these corrections very nearly balance one another 
in the case .of Venus, though they do not do so with the other 
planets, nor do they in the case of Venus with arguments other 
than <f>. 

I have compared most of Radau’s individual results with mine, 
and find a close agreement for nearly all those portions which he has 
computed. For the direct action this agreement was particularly 
useful, since my coefficients are obtained from formulae and methods 
radically different from his, though reducible to them by algebraioal 
processes. 

Two other differences should be noticed, namely, those in the 
terms with arguments — D + 4T — 3 V and l - D + 4T - 3 V. These 
arise from the fact that the indirect action is a little greater and of 
the opposite sign to the direct action; Radau only takes the latter 
into account. 

8. The differences B - N. These have been marked with the 
letters Ej, E a , E 8 (explained), and U (unexplained), whenever they 
exceed o 7, /o20. 

Those marked E : are due to the fact that Newcomb substitutes 
the variations of the elements in the value of V limited to the 
three terms 

V= w x + 2e sin l + |e 2 sin 2I, 

instead, of taking the complete value; he adds that “in nearly or 
quite all cases we may drop terms of the second order in e.” 

This third term of V in reality only produces a few secondaries 
with coefficients greater than o" , o20, and these arise in the form 
2D - 2I + <f> + (2T) = 2D + <j>; but there are many greater than 
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o"*oo3, tlie adopted limit below which coefficients were dropped. 
The latter arise also id the forms l + <j> - (2I) and 2D — l + + (2/). 

A much more extended and larger set of differences arises, 
however, from the neglect of the evection in the expression for V. 
The majority of the differences E x due to this cause will be found 
in the combinations— 

(f> + (2D — T ), 
l -H (f) + (2D — Z) = 2D 4* <f >, 

2D - l + <f> — (2D - l) = <f>, 

2 1 — 2D + <f> + (2D 1 ) = l + <f >, 

the last being of importance only when 2I - 2D + </> is of long 
period of the order of ten years, and the first quite rarely. . Hence 
the final coefficients of terms with arguments <£, 2D + will be 
chiefly affected, and a very few in those with arguments l + <j>, from 
this cause. All those differences which may be nearly or altogether 
explained by the neglect of this and other terms in V are marked 
E 1# I assume that the terms not present in Newcomb's list have 
not been computed. 

In order to examine the differences E 2 a closer inspection is 
necessary. On referring to Newcomb’s tables for the variations of 
the elements (p. 154 of his work), I find that under the argument 
2D he has only the terms containing (V — T) and 2(V — T) for 
Venus. There are, however, several other arguments with coeffi¬ 
cients of the same order' of magnitude as these, and the same fact 
appears to a smaller extent with the other planets. All differences 
explained by this cause have been marked E 2 . 

The small difference E g in the great long period term due to 
Venus scarcely needs explanation. Newcomb has included the 
portion (less than o"‘io) due to the mutual perturbations of Venus 
and the Earth, and he states that his method indicates a possible 
error of the order o"‘io. I exclude the former part, but the 
maximum error by my method should be less than o"'c>5. 

Of the differences U there are two in Mars, eight in Jupiter, 
and one in Saturn. The most important is that in Jupiter with 
argument J ; I have compared the several portions of any value 
for the primary with that of Eadau, and the results (after the 
corrections noted in No. 7 above have been made) agree for the 
indirect action, but I have been unable to obtain his result for the 
direct action from the algebraical formulae which he gives. The 
causes for the remaining differences I have not been kble to trace, 
and must leave the values of those coefficients an open question. 

9. In working out the terms produced by planetary action, the 
valuable memoir of M. Eadau has been available for comparison of 
results at almost every stage, and has materially assisted in the 
prevention of errors made during the course of a piece of work of 
great complexity, though not of great difficulty when once the 
theory had been put into final shape. The work of Professor 
Newcomb only appeared when I had finished the greater part of 
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the computations, and would in any case not have been available 
for detailed comparison without much labour, owing to the complete 
difference between his method and mine. He combines the direct 
and indirect actions at the earliest opportunity. I have kept 
them separate until each was fully completed. The comparison 
between his results and mine revealed one error in my work which 
affected the primaries due to the terms with argument <j> in E by 
about 5 per cent., and a few of the small secondaries from these 
terms about twenty per cent.; an error in the equations of varia¬ 
tions had almost no effect. These errors have, of .course, been 
corrected in the results given above. 

10. Like Professor Newcomb, I have also made an examination 
of the inequalities of the second order relative to the planetary 
masses, and have so far found nothing that could sensibly affect the 
motion of the Moon. I have found an additional portion to the 
term with argument l + 3T - 10Y of the order of o"-2, but the period 
of this term is so long that it would scarcely affect the observations 
within the degree of accuracy at present obtainable. The motion 
of the node of Venus also affects the term with argument 
l + 16T - 18V by a quantity of the same order of magnitude. 
The terms of the second order in the Sun’s motion, as given in 
Newcomb’s Tables of the Sun, 

8F'= - o"'265 cos (4Mj - 7T X + 3V9 - o"-02i sin^Mj -7^ +3V!) 

+ 376 cos (3 J x - 8MJ + 4T 1 !)-!- 5-18 sin (3J X - 8M X + 4^) 

(where the suffix denotes the mean anomalies of the planets instead 
of the mean motions), or ' 

+ o"*266 sin (x 0, igo^+3i°'8) (period 300 years) 

+ sin (o°*202# + 23i°*2) (period 1780 years), 

will produce inequalities with the same arguments due to indirect 
action having the approximate coefficients . 

- o"‘04 and — o"'9 

respectively. These results are sufficiently accurate for tabular 
purposes, but I shall give a more complete computation, with an 
examination of the whole effect of the terms of the second order. 

New Haven, Conn. U.S.A.: 

1907 November 26. 


Postscript .—Since this paper was sent in, I have computed the 
terms due to the motion of the ecliptic and have discovered a few 
new inequalities containing the arguments of the planets. The 
most important is one with a coefficient o"’2i due to Jupiter, and 
having a period of 280 years. 

1908 January 10. 
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